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Abstract 

With the HI detector at the ep collider HERA, D* meson production cross sections have 
been measured in deep inelastic scattering with four-momentum transfers Q 2 > 2 GeV 2 
and in photoproduction at energies around W 1P as 88 GeV and 194 GeV. Next-to-Leading 
Order QCD calculations are found to describe the differential cross sections within theoret- 
ical and experimental uncertainties. Using these calculations, the NLO gluon momentum 
distribution in the proton, x g g(x g ), has been extracted in the momentum fraction range 
7.5 • 10~ 4 < x g < 4 • 1CP 2 at average scales fi 2 = 25 to 50 GeV 2 . The gluon momen- 
tum fraction x g has been obtained from the measured kinematics of the scattered electron 
and the D* meson in the final state. The results compare well with the gluon distribution 
obtained from the analysis of scaling violations of the proton structure function F2. 
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1 Introduction 



The most precise determinations of the gluon momentum distribution in the proton have been 
obtained so far from the analysis of scaling violations of the proton structure function F 2 [II 0]. 
This method is however indirect in the sense that F 2 at low values of the Bjorken scaling variable 
xb actually probes the sea quark distributions which are related via the QCD evolution equations 
to the gluon distribution. The local behaviour of the structure function F 2 at a given value of xb 
depends on the gluon distribution x g g(x g ) in a rather wide range of values of the momentum 
fraction x g , and the analysis requires the assumption of a certain functional form of x g g(x g ), 
the parameters of which are then determined in a fit procedure. 

More direct determinations of the gluon density can be obtained by reconstruction of the 
kinematics of the interacting partons from the measurement of the hadronic final state in gluon- 
induced processes. Such direct measurements are complementary to the indirect analyses: al- 
though still limited in statistics, they are in principle more sensitive to local variations of the 
gluon distribution. They are subject to different systematic effects and provide an independent 
test of perturbative QCD. Direct gluon density determinations have previously been performed 
using events with J/^ mesons in the final state [0] and dijet events [Q]. 

The production of heavy quarks in electron-proton interactions proceeds, in QCD, almost 
exclusively via photon-gluon fusion, where a photon coupling to the incoming electron inter- 
acts with a gluon in the proton by forming a quark-anti-quark pair. This holds both for deep 
inelastic scattering (DIS) and for photoproduction where the exchanged photon is almost real. 
Differential charm photoproduction cross sections in the range of experimental acceptance were 
found to be reasonably well reproduced by such a description in Next to Leading Order (NLO) 
QCD [|5[], and measurements of the charm contribution to the deep inelastic proton structure 
function, F 2 C [p, also confirm this picture. 

Compared to the dijet case, smaller invariant masses of the partonic sub-process can be 
accessed in charm production. It is thus possible to extend the gluon density determination 
towards smaller momentum fractions x g . The charm data are statistically less powerful, but 
there is also much less background from quark-induced processes. 

Recently, Next-to-Leading-Order (NLO) QCD calculations of differential charm cross sec- 
tions at HERA have become available, both for DIS [J7|, |8p, and for the photoproduction re- 
gime [§]. With these calculations it becomes possible to determine x g g(x g ) directly. Here, 
measurements of differential D* production cross sections with the HI detector are presented 
in both regimes, and from these results the gluon density is determined in NLO. 

The theoretical framework and the reconstruction of the gluon momentum fraction will be 
described in more detail in the next section. In Section 3 the measurement of differential D* 
production cross sections will be presented separately for DIS and photoproduction. In each 
regime, the results will be compared with predictions based on the appropriate NLO QCD 
calculations. Section 4 will explain the method used to extract the gluon density from the 
measured cross sections, and will present the results. The paper concludes with a comparison 
to results obtained by HI from the analysis of the structure function F 2 . 
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2 Principle of the Analysis 



2.1 NLO QCD calculations 

Differential charm production cross sections in ep interactions have been calculated [ JTOj , [TTJ] in 
NLO QCD in the MS renormalization scheme. The computations were carried out using the 
subtraction method, which is based on the replacement of divergent soft and collinear terms in 
the squared matrix elements by generalized distributions. Charm quarks are treated as massive 
particles appearing only in the final state whereas the three light flavours (and gluons) are the 
only active partons in the initial state ("massive" scheme, "Three Flavour Number" scheme). In 
Leading Order, the production proceeds via the partonic sub-process jg — > cc. In NLO, there 
are in addition contributions to the cross section associated with diagrams where a real gluon is 
radiated from a gluon or a quark line, 7$ — > ccg, there are quark induced contributions via the 
processes ^q — > ccq and 'jq — > ccq, and there are virtual corrections to the cross section. 

For the photoproduction processes, resolved photon-proton interactions have to be taken 
into account in addition to the direct photon-proton scattering. While in the latter processes, 
the photon couples directly to a parton of the proton, the resolved photon processes involve 
the partonic structure of the photon itself. In NLO, the distinction of these processes is not 
unambiguous; here, a part of the photon parton density is contained in the direct photon calcu- 
lation. The remaining part that depends on the parton density functions of the photon is termed 
resolved. In the calculations of charm production performed in the 'massive' scheme the direct 
contribution is dominant in the phase space where at least one charmed meson is produced in 
the detector acceptance. Previous measurements of charm photoproduction [|]] have confirmed 
this. 

Depending on the kinematic regime under consideration, the QCD scales in the theoretical 
calculations are governed by different quantities. In the DIS case the factorization scale is set 
to be fip = ^/4m^ + Q 2 , where m c denotes the charm quark mass and Q 2 the virtuality of the 
exchanged photon, and the renormalization scale is set to be /jr = fxp- In the calculation of 7p 
cross sections, the factorization scale is given by the transverse momenta of the charm quarks: 
I~lf = \J^m 2 + A.p\, and the renormalization scale is chosen to be hr = Hf/Z . Variation of 
these choices will be one means of estimating the uncertainties of the perturbative calculations. 

Since the computations are differential in all relevant quantities, they provide access to the 
four-momenta of the outgoing partons. The evaluation of cross sections can be performed using 
Monte Carlo integrators. The calculations by Harris and Smith for DIS have been implemented 
in the HVQDIS program pSp , and those for the photoproduction case by Frixione et al. in the 
FMNR package [§. 

Charm quark hadronization into D* mesons is performed in the programs by convoluting 
the charm quark cross section with the Peterson fragmentation function [O]; the D* three- 
momentum vector is formed by scaling the quark three-momentum vector, and the D* energy 
is fixed such that the mass is m(D*) = 2.01 GeV. Evolution of the fragmentation function, 
which one expects to become important for transverse c quark momenta pj_ ^> m c , is not 
included because with present experimental cuts the region of interest is p± ~ m c . Cross 
sections for charm photoproduction at large transverse momenta p± 3> m c have been computed 
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in the"massless" scheme [ |T3| , Bp , where charm acts as an active flavour in the initial state proton 
and photon. 



Using the HVQDIS [|]] and FMNR ^\ programs, cross sections incorporating experimen- 
tal cuts can be calculated in NLO as a function of any final state variable depending on the 
D* and electron four-momenta. Thus the predictions of the NLO QCD calculations can be 
confronted directly with measured differential cross-sections for inclusive D* production in the 
experimentally accessible kinematic range. 

2.2 Reconstruction of the gluon momentum fraction 

The kinematics of the Leading Order 2^2 process — > cc is completely determined if the 
momenta of one incoming and one outgoing particle, here the photon and the charm (anti-) 
quark, are known. The momentum fraction x g in the infinite momentum frame is given by 

Xa = Ml±Ql with M* = #±4 and ls £^= <*-»tt* , (1) 
y ■ s z(l- z) p-q 2yE e 

where p\ c is the transverse momentum of the charm quark with respect to the photon-proton 
axis in the photon-proton centre-of-mass frame. The Lorentz invariant z can be evaluated in 
the laboratory frame from the charm quark's energy E and its momentum component p z ^\ The 
variables p, q, and p c denote the four-momenta of the incoming proton, exchanged photon, and 
outgoing (anti-)charm quark, respectively, so that the usual DIS variables read Q 2 = —q 2 and 
y = (2p ■ q)/s; ^fs is the ep centre of mass energy, and E e the incoming electron energy in 
HERA. In DIS, Eq. [1] can be written in the form x g = x,b(1 + M 2 /Q 2 ), which relates x g to the 
Bjorken scaling variable i B . 

In the presence of gluon radiation (or a gluon transverse momentum in the initial state), 
relation (JT]) holds only approximately. Also, in the experiment one measures not quarks but D* 
mesons that on average carry only a fraction of the energy of the primarily produced c quark. 
However, one can construct a hadronic observable x° BS that is well correlated with the true x g 
by replacing in Eq. [j] 

Pie -1.2^. , (E-p z ) c ^(E-p z ) D * . (2) 

The factor 1 .2 is introduced for convenience in the unfolding procedure described in Sect. ^; it 
is chosen such that x° BS ~ x g on average. 

2.3 Analysis procedure 

The analysis proceeds in two main steps. First, the D* cross section is measured in the region 
of experimental acceptance as a function of various kinematic variables and of x° BS . The data 
are corrected only for detector effects. As there is no extrapolation into unmeasured regions of 

'HI uses a right-handed coordinate system with the z axis defined by the incident proton beam, and the y axis 
pointing upward. 
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phase space nor are there any parton-hadron correlations involved, this step can be performed 
using the Leading Order Monte Carlo programs which describe the final state reasonably well. 
The resulting cross sections are then confronted with the NLO QCD predictions. 

In the second step the cross section as a function of x° BS is unfolded to the true variable 
x g and the gluon density x g g(x g , n 2 F ) is extracted. This step is carried out in NLO QCD in the 
framework of the Three Flavour Number Scheme. The programs HVQDIS for DIS and FMNR 
for photoproduction are used to calculate the cross sections in the region of acceptance as a 
function of x® BS as well as the correlations between the hadronic and partonic quantities, x° BS 
and x g , taking into account NLO effects and fragmentation. 



3 Cross Section Measurement 



3.1 Experimental conditions 

The data used here have been collected in the years 1994 to 1996 with the HI detector at the 
HERA collider, where 27.5 GeV positrons (henceforth generically termed "electrons") collided 
with 820 GeV protons, at a centre of mass energy of ^fs = 300 GeV. 

The HI detector and its trigger capabilities have been described in detail elsewhere []T5j]. 
Charged particles are measured by two cylindrical jet drift chambers [ jl6l , |T7j ], mounted con- 
centrically around the beam-line inside a homogeneous magnetic field of 1.15 Tesla, yield- 
ing particle charge and momentum from the track curvature in the polar angular range of 
20° < 9 < 160°, where 9 is measured using the central jet and two z drift chambers and is 
defined with respect to the incident proton beam direction. One double layer of cylindrical 
multi-wire proportional chambers (MWPC) iTTSp with pad readout for triggering purposes is 



positioned inside and another one in between the two jet chambers. The tracking detector is 
surrounded by a fine grained liquid argon calorimeter [jT9|], consisting of an electromagnetic 
section with lead absorbers and a hadronic section with steel absorbers. It covers polar angles 
between 4° and 154°. In the backward region (153° < 9 < 177.8°) it is complemented by a lead 
scintillator "Spaghetti" calorimeter (SpaCal) [EDfl which is optimized for the detection of the 



scattered electron in the DIS kinematic range under consideration here and provides Time-of- 
Flight functionality for trigger purposes. Installed in 1995, it consists of an electromagnetic and 



a more coarsely segmented hadronic section. A four-layer drift chamber (BDC) Q2l| ] mounted 
on its front improves the angular measurement of the scattered electron. 

The luminosity is determined from the rate of Bethe-Heitler ep — > ep7 bremsstrahlung 
events. The luminosity system consists of an electron detector and a photon detector, located 
33 m and 103 m from the interaction point in the electron beam direction, respectively. The 
system is also used to tag photoproduction events by detecting electrons scattered at small 
angles. In 1995 and 1996 a second electron tagging detector of the same type was added at 
44 m. It covers a much lower y-range. Therefore the two electron detectors cover different 
regimes of photon-proton centre-of-mass energies, W 1P , and thus are complementary. Time- 
of-flight (TOF) systems in the forward and backward directions are used to reject beam gas 
background. 
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3.2 Extraction of the D* signal 



The reconstruction of D* mesons in the DIS and photoproduction regime follows closely the 
method described in the previously published measurement of D* production It makes use 
of the D* tagging technique [Q], i.e. of the tight kinematical conditions in the decay chain 
£)*+ _> D°n+ ow -> (A^7r+)7r+ J] for which the overall branching fraction is 2.62% ||23|]. 
Particle tracks passing some quality cuts and fulfilling transverse momentum cuts depending on 
the kinematic range, are first combined in unlike- sign charged pairs. No particle identification 
is applied. Assuming one track to be a kaon and the other to be a pion, the invariant mass is cal- 
culated for all combinations and each possible mass hypothesis assignment. Those K~n + pairs 
with an invariant mass consistent with the D° mass within ±80 MeV are combined with a third 
track ("7rJ ou ,") which has to have charge of opposite sign to that of the kaon candidate and to 
which is assigned the pion mass hypothesis. D* + production is found as a distinct enhancement 
in the distributions of the mass difference AM = M(K~ir + 7rf low ) — M(K~~n + ) around the 
expected mass difference of 145.4 MeV [j23|]. The AM distributions for all D* candidates are 
shown in Fig. [I] for the DIS (a) and the two photoproduction samples (b,c). No enhancement is 
observed if the mass difference for the wrong charge combinations M(K~n~n + ) — M(K~n~) 
is used instead. 

The number of events, Nd* , is extracted from maximum likelihood fits to the distributions of 
the mass difference with a function taken to be a superposition of a Gaussian for the signal and 
a term N B (AM — m n ) a for the background. The position and width of the Gaussian are fixed 
to the values found from the total sample including both photoproduction and DIS events. The 
position is consistent with the expected value, and the width agrees with simulation results. The 
normalization of the signal, the background normalization Nb, and the shape exponent a are 
free parameters in the fit. Uncertainties from the parameterization of signal and background are 
accounted for in the systematic error. From Monte Carlo simulations the contribution to the sig- 
nal due to D° decays into channels other than K~n + (mostly D° — > K~K + and D° — > n~n + ) 
is obtained to be r = (3.5 ± 1.5)%. For these reflections a correction is made. 



3.3 Acceptance and efficiency determination 

A Monte Carlo simulation is used to determine the detector acceptance and the efficiency of the 
reconstruction and the selection cuts. Electroproduction events and direct and resolved photo- 



production events were generated in Leading Order with the AROMA 2.2 fl24| ] and PYTHIA 



5.7 programs J2~5|]. The hadronization step was performed according to the Lund string model. 



The generated events were then processed by the HI detector simulation program, and were 
subjected to the same reconstruction and analysis chain as the real data. The observed rapidity 
and transverse momentum distributions of D* mesons are acceptably reproduced. The changes 
found when varying the shape of the charm fragmentation function by using a Peterson function 
with e c ranging from 0.026 to 0.046 are accounted for in the systematic error. The dependence 
of the simulated efficiencies on other parameter choices made for the simulation (charm mass, 
parton density distributions) was studied, found to be small, and is neglected. 



2 The charge conjugate is always implied. 
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Contributions from bb production, with subsequent decays of b flavoured hadrons into D* 
mesons, have been calculated using the AROMA generator to be about 1.6% for DIS and even 
smaller for photoproduction. No subtraction is made; the quoted D* cross sections are inclu- 
sive. A systematic error accounts for the change of the efficiency due to a variation of the bb 
contribution within a factor of 5 of its predicted value, thus taking the preliminary result p6] ] 
into account, where an unexpectedly high bb cross section at HERA was found. 



3.4 Differential cross sections 

Differential cross sections as functions of kinematic variables are measured by dividing the data 
sample into bins and determining the numbers of D* mesons and the efficiencies in each bin 
using the Am fit and correction procedure as outlined above. The effects of limited detector 
resolution on the variables are small with respect to the chosen bin sizes and are corrected for 
using the Monte Carlo simulation. 

The cross sections are determined in the ranges specified in Tab. [j], where r] = — In tan(#/2) 



Table 1 : Kinematic ranges of the cross section measurements. 



DIS 


7P, ETAG44 


7P, ETAG33 


2 < Q 2 < 100 GeV 2 

0.05 < y < 0.7 
p±(D*) > 1.5 GeV 
\r}(D*)\ < 1.5 


Q 2 < 0.009 GeV 2 
0.02 < y < 0.32 
p±(D*) > 2 GeV 
\y(D*)\ < 1.5 


Q 2 < 0.01 GeV 2 
0.29 < y < 0.62 
p±(D*) > 2.5 GeV 
\y(D*)\ < 1.5 



denotes the pseudo-rapidity in the laboratory frame. In photoproduction, where the photon- 
proton centre-of-mass frame moves along the z axis, the rapidity in the laboratory frame, 
y — — | In -f^S is conveniently used. 

In the following, the DIS and photoproduction analyses are presented separately, since the 
scattered electron is measured in different detectors and since different methods are used to re- 
construct the event kinematics. Different cuts in the D* selection reflect different combinatorial 
background levels in the kinematic regimes under consideration. 



3.4.1 Deep Inelastic Scattering analysis 

The data have been taken in the years 1995-96 after the installation of the SpaCal and the 
BDC. They correspond to a total integrated luminosity of C = 9.7 pb _1 . The events have been 
triggered by an electromagnetic cluster in the SpaCal of at least 6.5 GeV energy, in coincidence 
with a charged track signal from the MWPC and central drift chamber trigger. 

The identification of electrons is similar to the procedure used in the inclusive structure 
function measurement [E7p. Scattered electrons are identified as clusters in the SpaCal with 



energy E e > > 8 GeV, with a cluster radius <3.5 cm consistent with electromagnetic energy 
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deposition, and with a cluster centre matched by a charged track candidate in the BDC within 
2.5 cm. The scattering angle is required to be 9 e i < 177°. 

The DIS kinematic variables are reconstructed using the "electron (e)" and the "eE" me- 



thod [28]. In both methods 



9 i 

Q 2 = 4E e E e , cos 2 (3) 



whereas 



y e = l- {E e ,/E e ) sin 2 (^/2) and y eS = S , (4) 

where the sum in the denominator runs over all (liquid Argon and SpaCal) calorimetric signals 
as well as central drift chamber tracks, avoiding double counting. The sum in the numerator 
is over the same objects except those forming the scattered electron candidate. The results 
obtained with both methods are found to be in very good agreement. The eS method needs 
smaller QED radiative corrections, and it is less sensitive to uncertainties in the electron energy. 
The absolute energy scale of the SpaCal was known to a precision of ±2% (at 27 GeV) at the 
beginning of its operation in 1995, and to ±1% in 1996. The eS method is therefore used for 
the 1995 data. For the 1996 data, the e method is used, which yields slightly improved average 
resolution in x° BS . 

For D* reconstruction, cuts of p± > 300 (120) MeV are imposed on the on the transverse 
momenta in the laboratory frame K and ir (itsimu) candidates. In order to suppress combinatorial 
background, the range of z (Eqs. is restricted to z > 0.2 for p±(D*) < 3 GeV. Variations 
found by applying or not applying this cut are accounted for in the systematic error. 

A total of N D * = 583 ± 35 events is observed; a correction is made for a small fraction r 



which is attributed to reflections (see Sect. pT2| ). The efficiency for electron and D* reconstruc- 



tion and selection is found to be e reco = 42% on average, ranging from 34% in the lowest to 
57% in the highest x° BS bin. The trigger efficiency is determined from data, by using indepen- 
dent trigger conditions. It is found to be e t ri g = 98.6% on average for the 1995 running period, 
and 83.3% for 1996 when tighter track conditions were applied. The dependence of e tr i g on the 
event kinematics is small and is taken from simulations. 

There is no indication of photoproduction background seen in the data; the electron energy 
and y% distributions are consistent with expectations for a pure sample of DIS events. From 
simulations of direct and resolved charm photoproduction, we expect less than 1 event within 
the cuts specified. No subtraction is therefore made and a 1% error attributed to this possible 
background source. 



Radiative corrections have been calculated using the program HECTOR [|29J. They amount 
to 5 ra d = 3% (7%) on average for the eS (e) method and are below 10% in all bins. The cross 
section for the given range in Q 2 , y, p± and r\ (Tab. [I]) is obtained as 

N n *(l - r) 

' 1 (5.75 ± 0.35 ±0.79)nb , (5) 



e trig e rec BR(D* — > Kirn) C(l + 8 rad ) 
where the first error is statistical, the second is systematic and is detailed in Tab. 0. One of 
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Table 2: Systematic errors of the cross section measurements. 





DIS 


IP 


Trigger efficiency 


+4.3, -3.3% 


6.5% 


Track reconstruction 


7.5% 


7J*-Fit 


+3.7, -4.4% 


6% 


Reflections 


1.5% 


6-production background 


1% 


7P background 


+0,-1% 




Radiative corrections 


2.9% 




Energy scale elm. SpaCal 


3.3% 




Energy scale Had SpaCal (±7%) 


< 1% 




Energy scale LAr (± 4%) 


< 1% 




Electron angle (±0.5 mrad) 


< 1% 




BDC hit finding 


2% 




Electron tagger acceptance 




6% 


Luminosity 


1.5% 


branching fraction 


4% 


z - cut 


+7.5, -6.6% 




Fragmentation 


+0.6%, -2.4% 


3% 


Total 


+14.1,-13.4% 


15% 



the largest contributions in both analyses is due to the uncertainty in the track reconstruction 
efficiency. Other important sources of systematic errors include uncertainties in the trigger 
efficiency and in the D* signal extraction. In DIS, uncertainties related to the modelling of the 
D* distributions, as reflected in the dependence on the z cut, are also sizable. All contributions 
show only slight variation with x° BS . 

In Fig. |2| the cross section results are presented as function of p±(D*), p* ± (D*), r](D*), and 
Q 2 . Compared to the data are NLO QCD calculations performed with the HVQDIS program 
version 1.1 B using a charm mass m c = 1.5 GeV. The parton density set CTEQ4F3 [EOT] is 



used, which is consistent with the Three Flavour Number Scheme as implemented in HVQDIS. 
The value for the QCD scale, A^ 5 ^ = 237 MeV, corresponds in the two-loop approximation 
the world average for the strong coupling constant, a s (M z ) = 0.119 ± 0.002 [|2~3"|]. For the 
c — > D* fragmentation fraction, the value 27% is used, consistent with the average value of 
(0.71 ± 0.05)% for the whole fragmentation and decay chain c — > D* + — > K~n + ix + [ ]3~T| ] and 



the value of 2.62% Q23Q used for the decay branching fraction. For the fragmentation parameter 



e c = 0.036 is used, which has been extracted []3"2|] from ARGUS e + e annihilation data using 



NLO QCD in a Fixed Order massive charm scheme. 

The NLO QCD predictions (displayed as histograms) agree with the measurements except 
for da/dr] for which the calculations tend to undershoot the data at low r/ and to overshoot 
at high i]. This disagreement can be accommodated within acceptable ranges of parameters 
such as the QCD scale and parton density functions. The shaded band in the figures represents 
the uncertainty of the theoretical prediction arising from different choices of the charm mass 
between 1.3 and 1.7 GeV. 
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3.4.2 Photoproduction analysis 



Data are analyzed independently for the case where the scattered electron is detected in the 
electron tagger at 33 m and in the tagger at 44 m. Henceforth, the respective data samples will 
be referred to as ETAG33 (recorded in 1994-96) and ETAG44 (1995-96). The total integrated 
luminosities of these samples are 10.7 ptT 1 for ETAG33 and 10.2 pb -1 for ETAG44. The 
two samples have no events in common. The kinematic ranges within which the differential 
distributions are measured are given in Tab. |T]. 

The events were triggered by a coincidence of the ETAG signals with track candidate signals 
obtained from the MWPC and central drift chamber trigger systems. Proton beam induced 
background is reduced by excluding events with energy flow only into the forward region of the 
detector. 

The acceptances of the electron taggers and their trigger efficiencies are accounted for 
by assigning weights to individual events. The acceptance depends on y, which is recon- 
structed here using y = 1 — E e >/E e . For the ETAG33 sample it is above 20% with an average 
value of about 60% in the specified range. This corresponds to a ^p centre-of-mass energy 
range of 162 GeV< W 1P < 234 GeV, with a mean of (W 7P ) = 194 GeV and an average 
(Q 2 ) ~ 10~ 3 GeV 2 . The acceptance of the ETAG44 peaks sharply around y = 0.09 and corre- 
sponds to an average of (W JP ) = 88 GeV. For ETAG44, an average value for the acceptance is 
used for all events. 

The K and n candidates used in the D* reconstruction are required to have transverse 
momenta p± > 0.5 GeV (ETAG33) or p ± > 0.35 GeV (ETAG44) respectively. A cut of 
p± > 0.15 GeV is imposed on the slow pion track. The method of equivalent event num- 
bers [E3p is used to determine the errors of the acceptance-weighted numbers of D* candidates. 



The weighted numbers of combinations found in the Am fits are 489 ± 92 (for ETAG33) and 
299 ±75 (for ETAG44). 

Using the efficiencies for the track component of the trigger condition, as well as the D* 
reconstruction and selection efficiencies as determined by Monte Carlo simulations, the cross 
sections for the specified ranges are found using Eq. ^. In the photoproduction case QED 
radiative corrections are negligible due to the requirement of no energy deposition in the photon 
detector of the luminosity system. Systematic errors (see Tab. 0) are similar to those of the DIS 
analysis, except for those related to the electron measurement. Here, uncertainties related to the 
electron tagger acceptance give a sizable contribution. The systematic errors are combined in 
quadrature to give ±15%. 

Following the Weizsacker-Williams Approximation (WWA) [ft4|, the electroproduction cross 
section a ep is converted into a photoproduction cross section according to a(jp — > D*X) = 
a(ep)/F, where F is the flux of photons emitted by the electron. For ETAG33 (ETAG44) the 
flux factor is F = 0.0128 (0.0838). 

The measured distributions da{^p — > D* ± X)/dy are shown for ETAG33 and for ETAG44 
in Fig. |]a and c, respectively, for the range specified in Tab. |T|. The distributions da lp /dp± are 
shown in Fig. 0b,d. In addition, for ETAG33 the double-differential distribution d 2 a 7P /dydp± 
is presented for three different ranges of transverse momentum in Fig. |j. 
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The superimposed histograms shown in Fig. |3| and in Fig. |] represent the absolute predic- 
tions of the NLO QCD calculation [Q, using, as in the DIS case, m c = 1.5 GeV, e c = 0.036, 
and 27% for the c — > D* fragmentation fraction. The calculations are shown for the parton 
densities MRST1 [ j3~5| ] (proton) in combination with GRV-HO [ |3~6| ] (photon). The histograms 
for the ETAG33 sample are averages of calculations made at three representative W 1P values, 
weighted by the photon flux integrated over the represented range. For ETAG44 the calculation 
is performed at a fixed W 1P = 88 GeV. Reasonable agreement is observed for both the shape 
and the absolute normalization between the measured single differential pj_ and y distributions 
and the NLO QCD calculation. 

Fig. [| reveals in more detail that the agreement is good in the low p± region where the 
"massive" QCD calculation is expected to describe most reliably the data, and where the bulk 
of the events is found. For comparison, also overlaid in Fig. |]is a calculation performed in the 



"massless" scheme QT4|]. This approach is not expected to provide a description of the data in 
the region pj_ ~ m c . Fig. |] shows that the shape of the measured rapidity distribution in the 
upper p± range of this analysis is less well described in both approaches. 



4 Determination of the Gluon Density 

The gluon density x g g(x g ) is extracted from the measured cross sections da j ' d{\ogx^ BS ) which 
are compared to the NLO QCD predictions in Fig. ^. 

The photoproduction data (Fig. |5|b-d) are analyzed separately for three different photon- 
proton energies (88, 185, and 223 GeV). Here, the rapidity range is restricted to — 1.5 < y < 0.5 
(for ETAG33) and < y < 1 (for ETAG44). From the NLO calculation, the resolved photon 
contribution in this restricted range is expected to be less than 5% (using the GRV-HO p6| ] or 



LAC1 Q37Q parton density set for the photon) in every bin and is accounted for in the systematic 
error. 

The determinations of the gluon density in the proton from the DIS and from the photopro- 
duction cross section follow the same principles. The calculation of the differential D* cross 
section in the experimentally accessible range a(x° BS ) can be written as 

a ( x OBS\ _ I a„ ,.2\ ,,2\ A /„OBS „ ,.2M , „ lJ>BS, 
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dx g [g{x g , /j, 2 ) • a(x g , /j, 2 )A(x° bs , x g , /j, 2 )] + a quark (x^ BS ) (6) 



where ji is the factorization scale, a is the partonic cross section, and A is the integration kernel 
which contains the effects of gluon radiation as well as fragmentation and also incorporates the 
limited kinematic range of the measured cross section. It receives the dominant contributions 
from regions where x g ~ x° BS . The quark-induced contribution to the cross section (in the 
range of acceptance) is denoted a quar k- The form of Eq. || is such that a determination of 
the gluon density in the proton can be made from the measured cross section by means of an 
unfolding procedure. 

Firstly, according to Eq. the quark-induced contribution as predicted in the QCD cal- 
culation is subtracted from the data bin by bin. It is about 1% in total in DIS. The relative 
contribution rises with x° BS , but remains below 10% in the highest DIS and below 20% in the 
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highest 7p bin. The correlation with the gluon density via QCD evolution can be neglected here, 
but the uncertainty arising from it is included in the error. 

The quark-subtracted cross sections are converted into cross section distributions as a func- 
tion of the true x g using the iterative unfolding procedure for binned distributions of Ref. p8[]. 
The effects of limited detector resolution on x° BS are small with respect to the bin size and 
have already been corrected for in the measured cross section. The correlation between x° BS 
and x g as obtained from the QCD calculation is shown for the DIS case as a two-dimensional 
histogram in Fig. ^. This correlation is used for the first unfolding iteration and re-weighted for 
each subsequent iteration using the result of the previous one. After 4 steps at most, convergence 
is reached, and the "smeared" distribution a(x° BS ) agrees with the measurement. 

The unfolded cross section distribution factorizes: 

<r{x Bii ) ~ g(x gti , f4) ■ cr(x g ,i, Hi) ; (7) 

but for every bin i in x g (and W) the gluon density is probed at a different factorization scale 
Hi which depends on the phase space region. The scale in the DIS measurement changes from 
17 GeV 2 in the first to 28 GeV 2 in the last bin in x g , the average being (fi 2 ) = 25 GeV 2 . 
The Q 2 dependence of the cross section is well reproduced (see Fig. ^d). In photoproduction 
the factorization scale varies between 30 and 140 GeV 2 , with an average value of 50 GeV 2 , as 
can be seen in Fig. [7| where g(x g ) as determined from the photoproduction data is plotted as a 
function of the scale \i associated to the x g and W range of each measurement. To determine 
g(x, (fi 2 )) at an average scale, all data points are scaled to an average (/i 2 ): 

/ / 2\\ _ g( x g,ii (At 2 )) 1 '^ 601 "^ g( x g,ii A j - 2 ) t ^ ieor ^ , \experiment ,a\ 
g\ x g,ii \H 1 1 i 2\ltienrv ' ( \thenrv ' ^ 9,i) • ^oj 

The change of g(x gt i,H 2 ) when varying the factorization scale from /i 2 to (/i 2 ) is taken from 
theory, assuming the same scale dependence as in CTEQ4 (DIS) or MRST1 (tp). This change 
is however small: in the DIS measurement g(x, Hi) deviates from g(x, (fi 2 )) by no more than 
9%. In photoproduction, the variation can be inferred from Fig. [7]. The effect is included in 
Eq. m and the dependence on the input parton density scale variation is accounted for in the 
systematic error. 

The resulting gluon distribution, x g g(x g ), is shown in Fig. [5]a for the DIS case at the av- 
erage scale of /i 2 = 25 GeV 2 . The covariance matrix for the statistical uncertainty has been 
calculated; the inner error bars in the figure represent the diagonal elements only. The correla- 
tion coefficients (defined as = cov^/ ^/coVjjCOV^) are 18%, 27%, 42% for the first, second 
and third pair of neighbouring bins, respectively. An overall normalization uncertainty arises 
because of the experimental systematic error of the D* cross section. It is increased by the addi- 
tional uncertainty from the c — > D* fragmentation fraction, and is added (in quadrature) to the 
statistical error to give the total experimental error shown as outer error bars. The downward ex- 
tension of the error bar in Fig. ||a represents the decrease of x g g(x g ) if the bb contribution to the 
D* cross section is 5 times the size predicted by the AROMA Monte Carlo. The corresponding 
result from the photoproduction analysis is shown in Fig. |8|b. 

Theoretical uncertainties affecting the gluon density extraction are determined by re-cal- 
culating with the QCD programs for different parameter choices the cross sections and the 
correlations needed in the unfolding procedure. 
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The dominant contribution in the DIS case arises from the uncertainty due to different 
choices of the charm quark mass between 1.3 and 1.7 GeV. In addition, the QCD factorization 
and renormalization scale /i was varied between 2m c and 2 J Ami + Q 2 m DIS, following [|8p. 

For the photoproduction results, the scale variation dominates the theoretical uncertainty 
over the whole range; a factor 2 is allowed in both directions. Here the charm quark mass 
uncertainty affects mainly the scale and is not added as a separate contribution. 

One of the larger uncertainties is due to our ignorance of the charm fragmentation function 
which has not yet been measured in ep scattering. It is assumed here that this can be taken from 
measurements in e + e~ annihilation. A fragmentation parameter e c range between 0.026 and 
0.046 is allowed in the calculation. Included in this range are the values extracted from ARGUS 
and OPAL data in a Fixed Flavour Number scheme [[32|]. The average uncertainty in x g g(x g ) 



due to this variation of the fragmentation parameter is 4% in DIS and 10% in photoproduction. 

The uncertainty in the strong coupling constant (world average) is taken into account by 
variation of the scale A^ 5 ^ between 213 and 263 MeV [p3|]. In order to assess the dependence 



on the input parton density parameterization, the MRS(A') [ p9| ] set is used. In photoproduction 

-4% 
1-0% 



an additional error of ,ty accounts for resolved contributions and the uncertainty in the photon 



structure function]]. 

The contributions to the theoretical error due to the discussed variations in the parameters of 
the NLO QCD calculations are added in quadrature and displayed as shaded histograms along 
the abscissa in Fig. || The errors at small x g are dominated by the statistical uncertainty. 

The gluon densities extracted from DIS and photoproduction cross sections are in very good 
agreement with each other, as can be seen in Fig. They can be compared with other results 
determined in NLO from other processes. For example, the result for x g g(x g ) determined by 
HI [[]]] from a QCD analysis of the scaling violations of F 2 is also overlaid in Fig. |9|. This result 
has been obtained from a fit to an inclusive measurement of the proton structure function F 2 
using the NLO QCD evolution equations and the boson-gluon fusion prescription for the charm 
treatment. Within errors it agrees well with the gluon distributions determined in the analyses 
presented here, which are performed in the framework of the Three Flavour Number Scheme 
for DIS as implemented in the HVQDIS program, and in the "massive charm" calculations for 
photoproduction, implemented in the FMNR program. This amounts to an important verifica- 
tion of our understanding of the application of QCD in low— x DIS and to a demonstration of 
the universality of the gluon distribution in the proton within the framework of such calculations 
at NLO. 



5 Conclusion 

Differential D* cross sections have been measured in deep inelastic ep scattering and in photo- 
production. The results have been presented as a function of various kinematic variables and in 
different kinematic ranges. They have been compared to NLO QCD calculations in the "mas- 
sive" Three Flavour Number Scheme. The calculations were performed using the programs 

3 In DIS, resolved contributions are expected to be negligible [J£(J. 
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HVQDIS and FMNR, which have implemented the Peterson form of the fragmentation func- 
tion as determined from e + e annihilation data to model the non-perturbative hadronization 
process. The overall description of the differential cross sections by the theoretical calculations 
has been found to be adequate. 

Using the explicit NLO QCD calculations, the gluon density in the proton, x g g(x g ) has been 
determined from the measured cross sections as a function of x® BS . The experimental observ- 
able x° BS is used as an estimator for the gluon momentum fraction x g using the kinematic 
information from the D* meson in addition to the scattered electron. The function x g g(x g ) has 
been extracted at an average scale fi 2 = 25 GeV 2 . The gluon densities obtained from electro- 
production and from photoproduction of charm agree well with each other and with that derived 
from the QCD analysis of the inclusive measurement of the proton structure function F 2 . 
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Figure 1: (a) Mass difference AM = M(K~tt + ttJ 1ow ) - M(K~n + ) distribution of D* candi- 
dates in DIS. The solid line represents the result of a fit as described in the text, which is used 
to extract the D* cross section. For the photoproduction samples ETAG33 (b) and ETAG44 (c) 
the unweighted AM distributions are shown. 
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Figure 2: Differential DIS cross sections in the kinematic range of experimental acceptance (see 
Tab. The HI data are shown as points with error bars (inner: statistical, outer: total); the NLO 
QCD prediction using the CTEQ4F3 parton distribution set and a charm mass m c = 1.5 GeV is 
shown as a histogram. The shaded band represents the variation of the theoretical cross sections 
due to different choices of m c between 1.3 and 1.7 GeV. (a) Transverse D* momentum in the 
laboratory frame, (b) Transverse D* momentum in the hadronic centre-of-mass frame, (c) D* 
pseudo-rapidity, (d) four-momentum transfer squared, Q 2 . 
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Figure 3: Differential photoproduction cross sections in the kinematic range of experimental 
acceptance (see Tab. [I]). The data are shown as solid dots with statistical error bars for the 
ETAG33 sample (a,b) and for the ETAG44 sample (c,d). A common systematic error of 15% 
is not shown. The histograms show the NLO QCD predictions based on the FMNR program, 
using the MRST1 parton density parameterization for the proton and the GRV-HO set for the 
photon. The width of the shaded band represents the variation of the theoretical cross sections 
due to different choices of m c between 1.3 and 1.7 GeV. 
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Figure 4: Double-differential cross section d 2 a JP /dydpj_ for three different ranges of transverse 
momentum at an average W 1P = 194 GeV, shown as data points with statistical error bars. 
A common systematic error of 15% is not shown. The histograms represent the NLO QCD 
predictions in the massive scheme based on the FMNR program (using the parton density sets 
MRST1 (proton) and GRV-HO (photon)); the shaded bands show their variations due to differ- 
ent choices of the charm quark mass between 1.3 and 1.7 GeV. The curve shows the result of 
a NLO QCD calculation in the massless scheme using the CTEQ4M parton density set for the 
proton and GRV-HO for the photon, and a charm fragmentation function extracted from e + e~ 
in the same framework. 
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Figure 5: (a) Differential DIS cross section as a function of x° BS in the kinematic range of 
experimental acceptance (see Tab. [I]). Figures (b-d) show the corresponding photoproduction 
cross sections, separately for W 1 p = 223, 185 and 88 GeV, respectively. The meaning of the 
data points, histograms and their errors is the same as in Fig. |2| (for DIS) and Fig. ^ (for "fp). 
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Figure 6: Correlation between the observable x° BS — calculated using Eqs. — and die true 
momentum fraction x g , as predicted in NLO QCD. The bin area corresponds to the cross section 
contribution from the x g interval to the observed x° BS range. The correlation is distorted due 
to the effects of gluon radiation and fragmentation. 
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Figure 7: Gluon density distribution, measured in photoproduction as a function of the scale fi 2 . 
Data are plotted as g(x g , fx 2 ) and not x g g(x g , ll 2 ), in order to show most clearly the evolution 
with scale. Each key symbol represents a different bin in x g with the average value given in the 
figure. The error bars represent the statistical errors only. The dotted lines represent the gluon 
density of the MRST1 parameterization, shown as a function of scale li 2 for the same set of x 9 
values. 
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Figure 8: Gluon density, determined from DIS (a) and from from photoproduction data (b), 
compared to the CTEQ4F3 parameterization at jj 2 = 25 GeV 2 (full line). The error bars re- 
present the statistical (inner) and total experimental (outer) error; in addition the effect of a 
subtraction of bb background (if 5 times higher than predicted in a LO Monte Carlo program) is 
displayed in (a), in (b) the bb contribution is negligible. The theoretical systematic error is shown 
as shaded band along the abscissa and is dominated in DIS (a) by the uncertainty associated with 
the charm quark mass, in photoproduction (b) it is dominated by the uncertainty associated with 
the renormalization and factorization scale. 
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Figure 9: Gluon densities obtained from the two D* analyses. The inner error bars represent 
the statistical and the outer the total error. The systematic error is a quadratic sum of all con- 
tributions, dominated by the theoretical uncertainty. Both results are compared to the result of 
the HI QCD analysis of the inclusive F 2 measurement jj^] at jj, 2 = 25 GeV 2 (light shaded band) 
and the CTEQ4F3 parameterization. 
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